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I.  INTRODUCTION 


High  current  racetrack  beam  Induction  accelerators  and  modified  beta¬ 
trons  are  a  subject  of  increasing  Interest  as  sources  of  high  power  elec¬ 
tron  beams  for  free  electron  lasers,  flash  radiography,  and  other  applica¬ 
tions.  The  racetrack  Induction  accelerator  geometry  Is  illustrated  in  Fig¬ 
ure  1.  The  beam  Is  injected  from  a  conventional  pulsed  diode  beam  genera¬ 
tor  into  the  drlfttube,  is  progressively  accelerated  as  it  repetitively 
passes  one  or  more  Induction  modules,  and  then  Is  extracted  from  the  accel¬ 
erator  for  Its  Intended  use.  Beam  extraction  may  even  be  unnecessary  for 
microwave  applications,  because  a  slow-wave  or  rlppled-magnetlc-field  cav¬ 
ity  can  be  Inserted  in  a  straight  section  of  the  drlfttube.1 

Most  beam  stability  studies  for  high  current  recirculating  devices 
have  dealt  with  negative  mass  and  resistive  wall  instabilities.2-5  How¬ 
ever,  experience  with  linear  Induction  accelerators  suggests  that  beam 
breakup  and  image  displacement  Instabilities  due  to  beam  interaction  with 
the  Induction  modules  and  other  discontinuities  In  the  drifttube  may  be 
significant.6*7  The  beam  breakup  Instability  arises  from  a  resonant 
coupling  between  beam  transverse  oscillations  and  m=l  electromagnetic  cav¬ 
ity  modes  localized  to  the  acceleration  gaps,8  while  the  image  displacement 
Instability  Is  caused  by  Interrupting  the  m»l  beam  Image  current  In  the 
drlfttube  wall.9  Clearly,  the  two  are  Interrelated;  in  some  contexts  the 

Image  displacement  Instability  can  be  viewed  as  the  low  frequency  limit  of  v 

% 

the  beam  breakup  instability. 

This  paper  extends  Instability  results  developed  for  high  current  lin¬ 
ear  Induction  accelerators  to  cyclic  devices.  The  primary  differences  are 
(1)  that  the  beam  passes  the  same  few  gaps  again  and  again,  and  (2)  that 
the  curved  beam  trajectory  couples  longitudinal  to  transverse  motion,  per¬ 
haps  leading  to  a  hybridization  of  the  negative  mass  and  beam  breakup 
modes.  Experience  with  Interaction  among  Instabilities  In  other  situations 
suggests  that  the  latter  Item  Is  the  less  Important  of  the  two.10 
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Figure  1.  Simplified  representation  of  cyclic  induction  accelerator  1 
racetrack  drift  tube,  acceleration  gaps,  and  injection  and 
extraction  ports. 


Therefore,  we  consider  only  the  former  aspect,  periodicity  of  the  inter¬ 
action,  and  omit  curvature  effects.  As  in  earlier  calculations,  we  assume 
that  the  beam  propagates  at  the  speed  of  light,  that  lateral  displacements 
of  the  beam  are  small  and  rigid,  and  that  the  gaps  are  narrow  compared  to 
the  beam  oscillation  wavelength.  The  consequences  of  relaxing  these 
approximations  for  the  image  displacement  instability  under  different  cir¬ 
cumstances  is  discussed  elsewhere.9  We  expect  the  three  approximations  to 
be  adequate  for  present  purposes. 

A  general  dispersion  relation  for  a  multiple  gap  racetrack  induction 
accelerator  is  derived  in  Section  II.  Because  the  resulting  expression,  a 
determinant.  Is  cumbersome,  we  then  specialize  to  the  case  of  identical, 
uniformly  spaced  gaps,  obtaining  approximate  analytical  growth  rate  form¬ 
ulas  for  the  beam  breakup  and  image  displacement  instabilities.  Section 
III  Illustrates  Instability  growth  for  a  1  kA  electron  beam  accelerated  to 
40  MeV  through  fifty  cycles  of  a  four  gap  device.  (These  parameters  are 
based  on  an  Induction  llnac  developed  by  the  National  Bureau  of  Standards 
and  now  being  modified  at  the  Naval  Research  Laboratory.11)  First,  growth 
rates  are  estimated  based  on  the  formulas  of  the  preceding  Section.  Then, 
more  precise  numerical  values  are  obtained  with  the  beam  transport  code 
BALTIC,  Including  the  effects  of  acceleration  and  initial  transients.  Our 
findings,  summarized  In  Section  IV,  are  encouraging:  The  Image  displace¬ 
ment  Instability  is  negligible,  and  the  beam  breakup  instability  is  manage¬ 
able  provided  that  the  mode  quality  factor  Q  Is  kept  low.  A  comparison 
with  the  negative  mass  and  resistive  wall  instability  growth  rates  also  is 
provided. 

II.  DERIVATION  OF  GROWTH  RATE  EXPRESSIONS 

In  the  long  wavelength,  paraxial  approximation  the  linearized  equation 
of  transverse  motion  for  the  beam  centroid  Is  given  by9>10 


(1) 


(If  *  I?)  t*  (|r  +  -  '“cl  s 

■  i  5  "  J  Fj  8(z '  2j' 5 

Here,  £  *  x  +  iy  Is  the  (complex)  transverse  displacement  of  the  beam 
centroid,  p  is  the  beam  density,  y  *s  the  beam  relativistic  energy,  a  is 
the  beam  radius,  R  is  the  drifttube  radius,  and  «c  is  the  cyclotron 
frequency  for  the  magnetic  guide  field.  Time  is  normalized  to  the  beam 
axial  velocity,  which  is  in  turn  assumed  to  be  equal  to  the  speed  of 

light.  The  last  term  on  the  left  side  of  Eq.  (1)  represents  the  combined 

transverse  forces  on  the  beam  from  Its  Image  charge  and  current  In  the 
smooth  drifttube  wall.  Discontinuities  in  the  wall  are  treated  as 
Impulsive  forces  appearing  on  the  right  side  of  Eq.  (1). 


We  begin  the  derivation  of  the  dispersion  matrix  by  assuming  the  beam 
parameters  to  be  constant  in  time  and  Fourier  transforming  Eq.  (1). 

(1u  -  (ia>  -  +1  l 
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As  mentioned  In  the  Introduction,  we  shall  take  account  of  time-varying  y 
In  Sec.  Ill  by  direct  numerical  Integration  of  Eq.  (1). 

Next,  Eq.  (2)  Is  solved  for  £,  treating  the  right  side  of  the  equation 
as  a  source  term.  This  can  be  done  using  a  Green's  function  or,  equiva¬ 
lently,  by  twice  Integrating  the  equation  and  Imposing  periodicity  at  z  = 
l,  the  accelerator  path  length.  The  result  is 
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Y(k+  -  M  5(z)  = 

(3) 

N 

(2sink+L/2)-1  l  Fj  ^(Zj)  exp[ik+(sign[zj  -  z]L/2  -  (z^.  -  z))] 

1  N 

-  (2sink_L/2)  1  j  Fj  ?(z..)  exp[ik_(sigri[zJ.  -  z]L/2  -  (z^  -  z))] 

where  N  is  the  number  of  acceleration  gaps  and  Sign(x)  is  the  sign  of  x. 
The  wavenumbers  k+  characterize  the  beam  cyclotron  and  spacecharge  waves, 
respectively,  as  determined  by  the  left  side  of  Eq.  (2),  l.e. ,  between 
gaps. 

k±  -  «  +  («c  ±  l <oc 2  -  2Pa2/YR2]1/2)/2Y  (4) 


For  typical  Induction  llnac  parameters,  k+  «  u  +  uic/y  and  k_  «  u. 
Although  it  may  not  be  immediately  apparent,  s(z)  as  expressed  in  Eq.  (3) 
is  everywhere  continuous;  its  first  derivative  is  discontinuous  at  zj. 


Evaluating  Eq.  (3)  at  each  of  the  gaps  yields  the  system  of  equations 
Y(k+  -  kj  5  (z£)  = 

(5) 

l  N 

(2s1nk+l/2)  1  l  Fj  c(*j)  exp[1k+(L/2  -  Mod[z.  -  z^L])] 

J  *"  1 

1  N 

-  (2s1nk_l/2)  1  l  Fj  e(Zj)  exp[1kjL/2  -  Hod[Zj  -  z£,L])] 


Here,  Mod  (x,L)  is  x  modulo  L,  defined  to  lie  between  0  and  L.  The  deter¬ 
minant  of  the  coefficient  matrix  of  Eq.  (5)  Is  the  desired  dispersion  rela¬ 
tion.  In  general.  It  must  be  solved  numerically  to  obtain  u.  Additional 
analytical  progress  can,  however,  be  made  when  the  gaps  are  spaced  uni¬ 
formly  and  all  have  the  same  response  function  F. 
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(6) 


Y(k+  -  k_)  5  UL/N)  * 

S(jl/N)  exp[ik+(L/2  -  Mod[j  -  i,N]L/N)] 

5(JL/N)  exp[1k_(L/2  -  Mod[j  -  £,N]L/N)] 

Since  Its  right  side  Is  a  discrete  convolution,  Eq.  (6)  Is  readily 
solved  by  performing  a  finite  Fourier  transform.  We  find 

y(k+  -  k_)  =  IF  {  (1  -  exp[2it1m/N  -  ik+  L/N])'1 

1  ,  W 

-  •(!  -  exp[2irim/N  -  ik_  L/N])  1  } 


(2s1nk  L/2)”1  F  ? 

j*0 


-  (2s1nk  L/2) ” A  F  Y 

j*0 


or,  more  simply, 

F  k  L/2  -  urn  k  L/2  -  irm 

1  *  '2y'(k+  V  )  i  ctnf - R - 1  "  ctnt - R - ]  1  * 

The  transform  index  m  ranges  between  0  and  N-l.  Still  another  useful 
representation,  reminiscent  of  results  for  nonrecirculating  devices, 
is8 


cos[  (k+  +  k_)L/2N  -  itm/N]  « 

(9) 

cos[(k+  -  k_)L/2N]  +  y(Tc"~Tc")  sin[  (k+  -  k_)L/2N]  . 
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For  F/wc  small  these  equation  are  approximately  satisfied  by  k+  *>  kn  or 
k_  -  kn,  with  kn  =  2irn/L.  Correspondingly, 

«  -  «n  =  kn  -  (<uc  ±  C^2  -  2Pa2/YR2]1/2)/2Y  .  (10) 

Instability  may  occur  when  u>  matches  a  resonant  frequency  of  F  or  when  the 
right  side  of  Eq.  (9)  exceeds  unity  away  from  a  resonance.  These  two  pos¬ 
sibilities  are  realized  in  the  beam  breakup  and  image  displacement  insta¬ 
bilities,  respectively. 


The  beam  breakup  instability  is  caused  by  TM^g  oscillatory  elec¬ 
tromagnetic  fields  localized  to  the  acceleration  gaps  and  adjacent  drift- 
tube  regions.  These  cavity  modes  are  conveniently  represented  by  damped 
harmonic  oscillator  equations  for  their  normalized  vector  potentials. 
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Q  is  the  cavity  mode  quality  factor,  while  u>0  is  its  frequency.  Coupling 
between  the  beam  and  a  cavity  mode  is  given  by  its  transverse  impedance 
Zj/Q.  The  transverse  impedance  is  basically  a  geometrical  factor,  which 
tends  to  scale  linearly  with  the  gap  width. 


Fourier  transforming  Eq.  (11)  in  time,  we  obtain  for  F  =  A/? 


F 


+  1wa>0/Q  -  u>0 
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The  beam  breakup  Instability  growth  rate  now  can  be  approximated  in  the 
usual  manner  by  expanding  Eq.  (8)  and  (12)  about  u»n,  defined  in  Eq. 
(10).  The  resulting  quadratic  In  oj  -  a»n  shows  instability  only  for  kn 
«  k_,  with  growth  rate 


Note  that  this  derivation  fails  for  sin[(k+  -  k_)L/2]  too  small;  see  Eq. 
(9).  More  precisely,  the  Instability  growth  rate  is  reduced  whenever 
| Nkn  -  a>c/y J  Is  less  than  the  absolute  value  of  the  right  side  of  Eq. 

(13)  and  finally  vanishes  when  Nkn  =  uc/y. 

The  Image  displacement  instability  arises  due  to  interruption  of  the 
beam  image  current  at  a  discontinuity  in  the  drifttube  wall,  such  as  occurs 
at  an  acceleration  gap.  The  analysis  in  Ref.  9  and  12  gives  a  force  coef¬ 
ficient 

F  *  K  (14) 

R*  e 

The  effective  gap  width  i  Is  given  by  the  physical  width  when  the  latter  is 
small  compared  to  the  drifttube  radius  R.  In  the  opposite  limit,  i  <  R. 
As  already  noted,  growth  occurs  whenever  the  right  side  of  Eq.  (9)  exceeds 
unity.  For  F/wc  small  this  happens  only  in  a  narrow  parameter  region 
centered  on  tan[(k+  -  k.)L/2]  =  F/y(k+  -  k. ) ,  for  which  the  growth  rate  is 
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r  »  NF/o>cL. 


(15) 


In  the  next  section  we  evaluate  these  two  growth  rates  for  possible 
recirculating  Induction  accelerator  parameters  and  then  generalize  our  num¬ 
erical  results  to  Include  Irregular  gap  spacing  and  beam  acceleration  by 
computer  solution  of  Eq.  (1). 

III.  NUMERICAL  EVALUATION  OF  GROWTH  RATES 

The  four  gap  linear  Induction  accelerator  at  the  Naval  Research  Labor¬ 
atory  Is  an  attractive  candidate  for  a  recirculating  device.  The  accelera¬ 
tion  gaps  are  In  pairs  with  the  members  separated  by  about  30  cm.  The 
distance  between  the  two  pairs  Is  arbitrary  within  reasonable  limits.  Here 
we  take  the  distance  to  be  200  cm,  giving  a  total  round-trip  path  length  of 
L  =  460  cm.  (The  specific  value  chosen  does  not  strongly  Influence  our 
conclusions.)  The  principle  m=l  gap  normal  mode  Is  at  880  MHz  with  a 
transverse  Impedance  of  15  ohms  and  a  Q  of  60.  Each  gap  provides  an  accel¬ 
eration  of  0.2  MeV,  resulting  In  an  energy  gain  for  the  electron  beam  of  40 
MeV  after  50  passes.  A  beam  current  of  1  kA  and  guide  field  of  2  kg  are 
assumed.  These  parameters  are  summarized  In  Table  1.  Note  that  y  In  the 
table  has  been  corrected  for  the  spacecharge  depression  In  the 
drlfttube  and  also  that  Zj/Q  has  been  expressed  in  dimensionless  form  by 
dividing  the  transverse  Impedance  by  30  ohms. 

Although  the  four  gaps  are  not  spaced  uniformly,  a  reasonable  upper 
bound  on  beam  breakup  Instability  growth  can  nonetheless  be  obtained  from 
Eq.  (13).  We  find  r  «  12.9*10“4  cm-1,  corresponding  to  about  30 
e-foldlngs  during  the  course  of  acceleration.  (We  assume  u0  =  u>n.) 
The  most  obvious  way  of  cutting  this  amplification  to  an  acceptable  level 
Is  by  reducing  Q.  Recent  work  at  Lawrence  Livermore  National  Laboratory  on 
short  pulse  induction  accelerators  suggests  that  a  quality  factor  as  low  as 


9 


TABLE  1.  NOMINAL  RACETRACK  INDUCTION  ACCELERATOR  PARAMETERS 


Path  Lengths 

L  =  460  cm 

Drlfttube  Radius 

R  =  7  cm 

Beam  Radius 

a  =  1  cm 

Guide  Field 

Bz  =  2  kg 

(*>c  =  1.173  cm-1) 

Beam  Current 

I  =  1  kA 

(v  «  0.0588) 

Beam  Energy 

U  «=  0.4  -  40  MeV 

(y  =  1.5  -  80) 

Number  of  Revolutions 

50 

Number  of  Gaps 

N  =  4 

Acceleration  per  Gap 

AU  *  0.2  MeV 

(Ay  *  0.4) 

Gap  Resonant  Frequency 

880  MHz 

Mode  Quality  Factor 

Q  =  60 

Gap  Transverse  Impedance 

15  ohms 

(ZJ/Q  *  0.5) 

Gap  Width 

X  -  5  cm 

10 


six  is  achievable.13  For  Q=6,  r  *  2.4«10-4  cm-1,  or  5.5  e-foldings.  As 
many  as  eight  e-fol dings  may  be  tolerable  provided  the  gaps  are  not  excited 
appreciably  at  their  resonant  frequency  prior  to  beam  injection. 

More  precise  growth  rates  can  be  obtained  by  numerically  solving  Eq. 
(5)  plus  Eq.  (12),  and  this  is  not  difficult  to  do.  Instead,  we  choose  to 
integrate  Eq.  (1)  plus  Eq.  (11)  directly  using  the  computer  program 
BALTIC.  The  code  was  exercised  extensively  in  support  of  the  RADLAC  radial 
pulseline  accelerator  program  and  so  is  well  tested.10  BALTIC  has  the 
advantage  over  a  dispersion  relation  solver  that  it  takes  account  of  beam 
acceleration  and  of  transients.  It  is,  of  course,  much  slower. 

Table  2  summarizes  thirteen  runs  of  the  BALTIC  code.  The  quantities 
N,  Q»  tu0»  and  y  were  varied,  but  with  NZ^/Q  held  fixed.  Some  calculations 
involved  an  accelerating  beam.  Cases  1-4  illustrate  the  effect  of  changing 
oj0  -  o>n.  Because  kn  =  0.014*n  cm-1,  the  pattern  of  growth  rate 
variation  with  a>0  in  these  four  cases  repeats  itself  with  the  same  0.014 
cm-1  periodicity  in  «0.  It  so  happens  that  sin[(k+  -  k_)L/2]  is  small 
for  y  »  80,  the  value  chosen  in  cases  1-4,  which  reduces  the  growth  rate 
some.  Setting  y  =  60  avoids  this  situation,  increasing  r  by  25%,  as 

illustrated  in  cases  5  and  6.  The  growth  rate  change  caused  by  dropping  Q 

from  60  to  10,  as  in  case  7,  is  consistent  with  Eq.  (13). 

The  seven  cases  just  discussed  were  for  a  single  gap.  Cases  8  and  9 

treat  two  gaps  evenly  spaced.  Growth  rates  are  comparable  to  the  single 
gap  runs,  as  expected.  The  enhancement  of  r  for  case  9  relative  to  case  8 
again  arises  from  adjusting  s1n[(k+  -  k_)L/2]. 

Finally,  cases  10-13  treat  four  gaps  spaced  as  described  at  the  begin¬ 
ning  of  this  Section.  Also,  the  beam  Is  accelerating  In  these  last  four 
runs.  It  Is  clear  from  the  resulting  growth  rates  that  the  gradual  accel¬ 
eration  which  occurred  has  little  effect  on  stability  except  to  average 


TABLE  2.  SUMMARY  OF  BEAM  BREAKUP  INSTABILITY  CALCULATIONS  WITH  BALTIC 


Case 

N 

Q 

“0, 

(cm-1 ) 

Y 

r 

(10-4  cm 

1 

1 

60 

0.18 

80 

8.3 

2 

1 

60 

0.1732 

80 

3.3 

3 

1 

60 

0.17757 

80 

8.9 

4 

1 

60 

0.1787 

80 

8.8 

5 

1 

60 

0.18 

60 

10.0 

6 

1 

60 

0.17757 

60 

12.3 

7 

1 

10 

0.18 

80 

3.4 

8 

2 

60 

0.18 

80 

9.6 

9 

2 

60 

0.18 

1.5 

12.0 

10 

4 

60 

0.18 

1.5-80 

9.8 

11 

4 

10 

0.18 

1.5-80 

2.6 

12 

4 

6 

0.18 

1.5-80 

1.2 

13 

4 

6 

0.17757 

1.5-80 

1.6 

12 


over  «0  -  wn  and  sin[(k+  -  k_)L/2],  In  case  13,  w0  was  decreased 
slightly  to  maximize  growth  at  large  y  Indeed,  r  Increased  by  one-third 
relative  to  case  12. 

Figure  2  depicts  the  maximum  transverse  displacement  of  the  beam  as  a 
function  of  time  for  case  13.  It  and  the  other  runs  were  initiated  with  a 
uniform  transverse  offset  of  unit  magnitude  for  the  beam.  The  cavity  modes 
were  Initially  excited  at  an  amplitude  consistent  with  injection  of  the 
offset  beam  with  zero  risetime,  a  worst  case  assumption.  The  displacement 
is  seen  to  grow  by  a  factor  of  about  3.6. 

The  Image  displacement  instability  turns  out  to  be  insignificant. 
From  Eq.  (15)  Its  peak  growth  rate  Is  0.9«10-4  cm-1.  Since  this  occurs 
only  over  a  narrow  range  of  parameters,  however,  we  should  expect  a  much 
slower  average  growth  during  acceleration.  Indeed,  a  four  gap  BALTIC  run 
gave  no  coherent  growth  whatsoever,  and  transient  fluctuations  amounted  to 
less  than  15%  of  the  Initial  displacement. 

IV.  CONCLUSIONS 

In  this  report  we  have  developed  a  simple  theory  of  beam  breakup  and 
Image  displacement  Instabilities  In  cyclic  Induction  accelerators  and 
applied  It  to  obtain  growth  rate  estimates  for  a  possible  device.  The 
theoretical  model  takes  recirculation  into  account  primarily  by  enforcing 
periodicity  on  the  unstable  transverse  models.  Drlfttube  curvature  effects 
were  not  considered.  We  find  for  the  1  kA,  40  MeV  accelerator  that  neither 
the  beam  breakup  nor  the  Image  displacement  Instability  should  be  serious. 
Specifically,  with  Q»6  and  parameters  otherwise  as  in  Table  1  the  beam 
breakup  growth  rate  Is  r  *  1 .6- 10“**  cm-1;  the  image  displacement  growth  is 
negligible. 


For  comparison,  growth  rates  of  the  negative  mass  and  resistive  wall 
(for  stainless  steel)  Instabilities  for  comparable  parameters  are  of  order 
5»10-4  cm-1  and  0.2* 10-4  cm-1,  respectively.3-5  Thus,  the  negative  mass 
Instability  Is  somewhat  more  serious  than  the  gap-induced  instabilities  in 
this  Instance. 

The  accelerator  path  length  L  *  460  cm  chosen  for  the  numerical  calcu¬ 
lations  Is  perhaps  small  for  practical  devices.  Increasing  L  by  lengthen¬ 
ing  the  straight  sections  of  the  racetrack  geometry  would  have  little 
Impact  on  any  but  the  resistive  wall  Instability  Integrated  growth:  Growth 
rates  of  the  gap-induced  Instabilities  decrease  roughly  as  L-1,  so  that 
growth  per  pass  remains  fixed,  while  the  negative  mass  Instability  grows 
only  in  the  curved  sections  of  the  drift  tube.  In  fact,  the  negative  mass 
mode  may  decay  In  straight  sections.  Growth  per  cycle  of  the  resistive 
wall  Instability,  of  course.  Increases  linearly  with  L.  Its  growth  rate  is 
so  slow,  however,  that  the  path  length  could  be  Increased  by  an  order  of 
magnitude  or  more  before  the  resistive  wall  Instability  became  competitive 
with  the  other  modes. 

The  Important  stabilizing  Influences  of  electron  energy  spread  and 
transverse  oscillations  have  been  omitted  from  our  calculations.  Recent 
estimates  for  betatrons  with  comparable  parameters  suggest  that  a  105 
spread  In  the  Injected  beam  electron  energies  Is  sufficient  to  suppress  the 
negative  mass  instability  provided  that  the  mean  Injection  energy  is 
greater  than  about  3  MeV.2-5  Additionally,  research  In  support  of  the  ATA 
accelerator  program  predicts  that  employing  occasional  nonlinear  focusing 
elements  to  spread  the  electron  transverse  oscillation  frequencies  can 
limit  beam  breakup  Instability  growth.14  We  remark  also  that  varying  the 
resonant  frequencies  n>o  from  gap  to  gap  can  cut  the  beam  breakup  growth 
rate  by  N-1 ,  provided  the  relative  shifts  are  greater  than  Q-1.15 


Extensions  to  the  present  research  naturally  fall  Into  two  categories, 
adding  drifttube  curvature  and  adding  beam  temperature.  The  former  is  rel¬ 
atively  easily  accomplished  by  inserting  the  transport  equations  of  Ref.  4 
or  16  into  BALTIC.  This  also  would  allow  us  to  determine  the  influence  of 
straight  drifttube  sections  on  the  negative  mass  instability.  Treating 
beam  temperature  beyond  what  has  already  been  done  probably  requires  numer¬ 
ical  simulation. 
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